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作用 下 质量 和 半径 的 演化 时 标 及 其 对 目 转 
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摘要 : 研究 了 带 有 引力 能 源 和 抛射 物质 的 引力 收缩 星 如 金牛 T 星 处 于 慢 引 力 收缩 阶段 
后 ， 在 抛射 物质 和 引力 收缩 联合 作用 下 对 质量 和 半径 的 改变 及 其 对 自转 角速度 变化 产生 的 影 
响 。 给 出 了 质量 和 半径 随时 间 演 变 的 联 立 微分 方程 组 及 其 解 。 利 用 解 计算 了 金牛 了 T 星 的 质 
量 和 半径 的 演化 时 标 以 及 对 自转 角速度 变化 产生 的 影响 。 给 出 数值 结果 ， 并 讨论 了 理论 和 数 
值 结果 。 
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恒星 形成 后 首先 进入 星际 云 快速 引力 收缩 阶段 ， 当 星际 云 的 内 部 压力 渐渐 增 大 时 处 于 准 流体 静 力 
平衡 ， 星 际 云 由 快 引力 收缩 阶段 进入 慢 引 力 收 缩 阶段 。 处 于 慢 引 力 收缩 阶段 的 恒星 发 生 不 规则 光度 并 
抛射 大 量 物质 ， 其 中 金牛 T 星 就 是 处 于 慢 引 力 收缩 阶段 抛射 物质 的 人 恒星。 在 快 引力 收缩 阶段 以 引力 
收缩 为 主导 ， 质 量 损失 可 不 考虑 , 因此 自转 加 速 加 快 。 当 进入 慢 引 力 收 缩 阶段 因 出 现 抛射 物质 时 质量 
损失 为 主导 而 引力 收缩 次 之 ， 因 此 自转 角速度 开始 减 慢 … 。 

对 于 金牛 T 型 星 抛射 物质 损失 质量 可 参见 文 [2-10]。 文 [2] 给 出 了 较 详细 的 研究 ， 但 星体 在 引力 
收缩 和 抛射 物质 作用 下 ， 质 量 损失 和 半径 的 演化 时 标 及 其 对 自转 角速度 变化 的 影响 并 没有 给 出 ,特别 
由 前 两 种 因素 作用 下 ， 质 量 损失 和 半径 收缩 对 金牛 T 型 星 的 自转 角速度 的 影响 没有 论述 。 

当 星 体 进行 引力 收缩 时 半径 缩小 ， 这 使 自转 角速度 加 快 ， 而 当 质 量 流 失 时 ， 自 转角 速度 变 慢 。 金 
牛 工 星 的 自转 角速度 变化 是 在 引力 收缩 半径 缩小 和 抛射 物质 使 质量 减少 这 两 种 因素 联合 作用 下 产生 
的 。 文 [11-16] 讨 论 了 有 关 金 牛 T 星 的 自转 变化 ， 但 较 少 同 半径 的 变化 相 联系 。 本 文 在 上 述 几 方 面 做 
了 进一步 研究 。 金 牛 T 星 是 每 个 恒星 诞生 后 早期 演化 必 经 的 过 程 ， 对 于 本 文 研 究 的 典型 金牛 T ER 
有 其 它 金牛 类 型 星 的 普遍 性 质 和 意义 。 


1 处 于 慢 引 力 收 缩 阶 段 金 牛 T 型 星 的 物理 模型 的 基本 假定 


(1) 金 牛 T 星 的 能 源 只 靠 引 力 收 缩 产生 的 引力 能 源 ， 尚 无 核能 源 ， 是 一 颗 带 有 引力 能 源 模型 的 
星 ， 故 星 的 质量 辐射 可 以 略 去 。 所 以 有 
1-B=P,/P=0, i B=1 
其 中 ，P 为 辐射 压力 ; PARA. 
(2) 星 的 质量 损失 来 自 某 种 物理 机 制 产 生 的 物质 抛射 ， 主 要 抛射 氢 粒 子 。 假 定 抛 射 物 质 的 速度 ， 
视 为 常量， 而 抛射 的 氧 粒子 的 密度 也 视 为 常量 。 
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(3) 星 的 半径 改变 主要 来 自 引 力 收缩 产生 的 变化 ， 其 次 也 与 抛射 物质 的 质量 动能 有 关 。 但 金牛 T 
星 抛射 物质 的 最 外 面包 层 是 膨胀 的 ， 而 内 层 半径 R 因 引 力 收缩 而 缩小 。 

(4) 金 牛 T 星 的 辐射 能 量 ( 热 光度 ) 主要 靠 自 身 引 力 收缩 释放 的 引力 收缩 能 和 抛射 物质 带 走 的 动能 之 差 。 

(5) 根 据 文 [17] 收 缩 星 在 引力 收缩 阶段 垂直 于 H-R 图 演化 ( Hayashi 轨迹 ) ， 即 演化 路 线 近 似 垂直 于 
H-R 图 上 的 横 坐 标 ( 有 效 温度 7.) ， 故 在 收缩 时 可 把 星 的 表面 温度 了 .近似 视 为 常量 ， 也 不 随时 间 变 化 。 

(6) 由 于 收缩 星 在 引力 收缩 阶段 只 靠 对 流传 送 能 量 ( 对 流星 ) ， 对 于 对 流星 其 多 方 指 数 取 n=1.5， 
y=5/3, 


2 决定 金牛 T 星 的 质量 和 半径 的 演化 方程 


处 于 慢 引 力 收 缩 抛射 物质 的 金牛 卫星 由 于 抛射 物质 损失 质量 ， 质 量 损失 率 可 由 下 式 确定 中 |. 

OT a 4aR?N,Vymy, (1) 
dt 

其 中 ，R 为 半径 ; Vy 为 在 表面 处 物质 的 抛射 速度 ， 它 可 由 观测 的 轮廓 给 出 ; Nu 为 氧 粒子 的 数 密度 ，; 
ma 为 氧 原子 的 质量 。 按 前 2 KRE, Vo, Ny 和 mi 为 常量 ， 而 m 和 RR NINA, Em, RI 
个 演化 方程 式 。 再 确定 m、R 的 第 2 个 演化 方程 式 。 根 据 星 的 构造 理论 "| ， 星 的 势能 E Sk 


4 
a 
E' =-B : 
2 y= 
3 m 
其 , Q=- Go 
ToO T a R 
oat E - 4/3 - 4/3 3 ; 
N 星 的 引力 收缩 能 : Be =- b=? pa=-* a are 
Ša y-1 y-1 5-n R 
Q 根据 第 1 节 第 (1) 和 第 (6) 条 假定 , B=1, n=1.5, y=5/3 
N 3 gm 
2 pv ASG (2) 


因为 星 收缩 时 只 能 引起 半径 的 改变 ， 不 影响 质量 改变 ， 所 以 收缩 能 随时 间 的 变 率 为 


dE 2 
c 3 ç” dR (3) 


© Ex = 二 mV,’. (4) 
所 以 带 走 的 动能 的 变 率 为 (5) 


按 第 1 节 第 (4) 条 假定 ,收缩 星 ( 金 牛 T 星 ) 的 辐射 能 量 Er 是 引力 收缩 能 Eo 和 抛射 物质 带 走 的 动能 Ek 


dE, dE, dE, o dE, 
HAZ, BR E,=E,-E,, Alt roor i 但 辐射 能 变 率 等 于 光度 了 ， 故 光度 了 等 
于 引力 收缩 能 减 去 抛射 物质 带 走 的 动能 : 
dE, dE | P 
7 ae o 
其 中 ， 热 光度 工 可 写成 : 
L=4rR oTt. (1) 


其 中 ，o 为 Stefan 常数 ; T, 为 星 表面 有 效 温度 。 
将 (3)、(5) 和 (7) 式 代入 (6) 式 中 ， 即 得 确定 m、R 的 第 2 个 演化 方程 式 : 


3 mdR 1 dm 
= = 4TR2o7 t. 8 
7 RR dt 2 dt tes (8) 


Ve 


ChinaXiv 合 作 期 刊 


DA 天 文 研究 与 技术 14% 


根据 第 1 节 第 (5) 条 假定 ， 收 缩 星 在 H-R ALE AIRE ee ar (SE EB A bs A Bit 7,)， 故 在 (8) 式 
右 端的 了 区 可 近似 视 为 常量 ， 又 根据 第 (2) 条 假设 ， 抛 射 物质 速度 了 也 可 视 为 常量 。 所 以 (1) 式 和 (8) 
式 组 成 的 质量 m、 半 径 R 随时 间 演 变 的 微分 方程 组 是 

dm 


qr 一 一 4 万 民 NTYomn， 
3 Cr2dR 1 dm a n4 
7 R a V 7 47R oT... 

再 将 第 1 式 代 入 第 2 式 后 得 金牛 T 星 在 引力 收缩 和 抛射 物质 作用 下 质量 和 半径 随时 间 的 演化 方程 组 : 
I 4aR?NymyV, . (9) 
dt 

dR 147R* 3 r 
dt 3Gm* (Nym Vy + 20T,") . (10) 


3 质量 和 半径 的 演化 时 标 ( 质 量 和 半径 同时 间 的 演化 关系 ) 


演化 方程 组 (9) ~(10) 是 可 积 的 ， 可 推出 质量 和 半径 的 演化 时 标 ， 要 求 得 质量 m 和 半径 R 随时 间 
变化 规律 和 数值 需要 进行 积分 求解 。 给 出 用 可 积 求解 法 所 得 到 的 质量 和 半径 的 演化 时 标 。 
由 (9) 式 和 (10) 式 消去 时 间 d 后 ， 首 先 可 得 质量 和 半径 在 演化 过 程 的 关系 式 : 


(n) 7 ( 207.° ， 
= +V, : 
R) dm 66 \Nam,V, 


变数 分 离 后 积分 : 
R dR _ ‘i dm 
Ro R i mo m?’ 
结果 有 R=m/(K, + K,m). (11) 
或 m= K,R/(1 - K,R). (12) 
7 20T* 
其 中 ， 6 = 本 | |= A. 13 
' 66\N,ym,V,  ° a 
1 K uy 
K,= 一 一 一 = 常数 . (14) 


Ry mo 
(11) ~(12) 式 是 质量 和 半径 演化 过 程 的 关系 式 。 首 先 求解 质量 的 演化 时 标 ， 将 (11) 式 的 RR 代入 (1) 
式 右 端 ， 有 


dm 4 m x y K m ? 
z- ir = 
dt K, + Km) Wto (K, + K,m 


K, + Km)? 天 + 2K Km + K, m? 
- K,dt = dm = dm. 
m m 
7 1 1 
积分 后 可 得 : - Ks(t -t)) =- Ki 二 - +) + 2K Kiln (=) + K,°(m — mo) 
m Mo Mo 
JM 1) =O 时 质量 的 演化 时 标 式 : 
K? Mo 
K,t = (m, — m)| K? + + 2K,K,In | 一 |. (15) 
mmo m 
其 中 K, = 4nNym Vo = 常数 (16) 


再 求解 半径 的 演化 时 标 ， 积 分 (10) 式 ， 令 : 
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ar = (NamaJ + 20T,*) = 和 常数. (17) 
dR Rt 
iu J 
上 有 dt t m? 


再 将 (12) 式 的 mm 代入 上 式 后 有 
dR (1 RR)R 
dt K, 2R? 


=K, (1 - K,R)’° R. 


积分 得 


t R dR 
Kf a =f G- KRR ; 


右 端 积分 式 可 用 分 解 有 理 分 式 为 部 分 分 式 的 积分 法 积分 ， 则 有 


R dR R2KdR re Ky dR R 2K,dR 
Kst =| 5 | | 2 + | _ ý 
R R R R Ro (1 - K,R) Ro | — K,R 


所 以 半径 演变 的 时 标 为 


1 1 K, (R, -R 1-K,R,)R 

Kt=(—-—)+ 2 (Ro R) pae o a (18) 
R, R’ (1-K,R)(1-K,R,) (1 - K,R)R, 

其 中 ， K, ]Kj/K (19) 


二 4 质量 损失 率 和 半径 收缩 率 对 自转 角速度 变化 的 影响 


= 当 星 体 引力 收缩 时 半径 缩小 ， 自 转角 速度 加 快 ， 而 当 质 量 流失 时 自转 角速度 变 慢 。 金 牛 T 星 的 
~N 自转 角速度 变化 是 在 引力 收缩 半径 缩小 和 抛射 物质 使 质量 减少 两 种 因素 联合 作用 下 产生 的 。 

= 当 星 体 抛射 物质 时 除 损失 质量 外 ， 还 要 损失 角 动 量 ， 损 失 的 角 动 量 是 由 抛射 物质 带 走 的 动量 。 所 
N 以 金牛 T 星 抛射 物质 时 损失 的 角 动 量 应 该 等 于 被 抛射 物质 带 走 的 角 动 量 。 设 金牛 了 星 的 自转 角速度 
> Ho, PEERK KR, TKR 是 抛射 物质 dm 的 回转 半径 。 如 果 抛 射 物质 是 各 向 同性 ( 球 对 称 )， 按 
角 动 量 守恒 有 020 


d( KR’ mo) = 天 Rodm . (20-1) 
= 写成 角 动 量 和 质量 损失 率 的 形式 : 
= d(KR’mo) dm 
© J = KR w EP (20-2) 


微分 后 有 
Ks-Kı ldm 2dR 1do_ 
K; m d Rd øw dt 


其 中 ，Ks 对 于 不 同 多 方 模型 星 的 指数 有 不 同 值 ， 当 恒星 进入 主 序 星 时 Ks 值 很 小 , 一 般 取 Ks=0.1， 
但 在 主 序 前 阶段 ， 如 金牛 了 星 的 K, 值 比 1 较 大 。 根 据 文 [20] 对 于 多 方 指数 n=1.5( 本文 所 取 的 值 为 
第 2 节 第 (6) 条 ) 所 对 应 的 K, E, ERTA Ks=1/5=0.2。 由 于 抛射 物质 各 向 同性 ， 所 以 取 K,=2/3. 
故 由 上 式 可 得 角速度 相对 变化 率 为 


. 7 R 
igen 20. (21-1) 
3m R 


如 果 考 虑 :=0 时 目前 的 角速度 相对 变化 率 ， 可 写成 : 

7 m(0) _ , R(0) 
3 m(0) R(O) ` 
w(0) 、m(0) 和 及 (0) 皆 为 目前 (:=0) 的 角速度 、 质 量 和 半径 ， 而 四 (0)、 交 (0) 和 尽 (0) 是 目前 (t=0) 
的 角速度 变化 率 、 质 量 损失 率 和 半径 收缩 率 。 如 果 考 虑 角速度 长 期 变化 ， 积 分 (21-1) 式 有 


ò (0)/w(0) = (21-2) 
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d 7 d 


—(] St ee 
ae) 3 dt 


d 
Inm) -2 二 (1 
di a ll p OR 


积分 dlnw = f dlnm 一 af dlnR 
wo 0 


mo 


7/3 2 7/3 R 2 
Im 
Wo mo Ry mo R 
m 1/3 R 2 
<. W/Wo = (= S) ; 
Mo R 


o =o, (m/m,)” (R,/R)?. (22) 
EF, w, m FRAN t=0 时 的 初始 值 。m 和 RR 随时 间 ; 变化， 由 (9)~(10) 式 确定 或 由 (15) 和 (18) 
式 确定 ， 故 角速度 w 也 是 时 间 i 的 函数 ， 随 时 间 变 化 。 如 果 将 (11) 式 的 RR 代入 (22) 式 后 可 得 角 速 w 


随 星 的 质量 变化 而 变化 ， 即 | 
m \3 K, + Km 2 
sof) EEE) » 
Mo kı + K,my 


5 理论 结果 对 金牛 T 星 (T Tauri) 长 期 演变 的 数值 计算 


利用 前 节 所 得 的 理论 结果 对 金牛 T 星 的 质量 、 半 径 和 自转 角速度 随时 间 的 演变 做 一 数值 计算 。 
先 用 文 [2] 给 出 金牛 T 星 的 各 项 物理 参数 . 

mo = 0.60me = 1. 193 4 x 10°g, Ry = 4.56R6, = 3.1737 x 10"'cm , 

T.(K)=4 100°K, N,=1.15 x10" 个 原子 /cem?，V=225 km/s 
按 物理 量 数 据 

M,,=1.6735x 10g, 7=5. 6695 x 10 erg/cm’ degts, G=6.67 x 10° (cgs 制 ) 
将 这 些 数据 代入 (13) 、(14) 、(16) 、(17) 式 得 

K =1.285 80x10” (cgs 制 ), K,=-1.088 75x 10° (cgs 制 ), K,=5.4414x 10° (cgs fil) , 

K,=7. 090 70 x 10" (cgs ți), K,=4. 175 60 x 10 (cgs tl) 

根据 以 上 数据 和 对 常数 天 计算 的 数值 ， 可 按理 论 式 子 推算 金牛 T 星 在 不 同 阶段 演化 时 标 对 应 的 
质量 、 半 径 和 自转 角速度 长 期 演变 的 数值 。 

首先 给 出 5 个 阶段 质量 损失 的 演变 数值 ， 即 m=0. 999 999 999m, 0. 999 999 90m。、0. 999 999 900m,、 
0. 999 999 000m。 和 0.999 990 000m。，(mo 是 现在 的 质量 ) 以 及 K, 和 ,的 数值 代入 (11) 式 可 得 5 个 演化 
阶段 的 半径 演化 数值 。 再 将 5 个 演化 阶段 的 质量 和 半径 的 演变 值 以 及 初始 值 m。 和 Ro 代入 (22) 式 得 5 
个 演化 阶段 的 自转 角速度 w 的 演变 值 。 最 后 再 将 初始 质量 m。 和 5 个 演化 阶段 的 质量 损失 数值 以 及 
Ki,、K, 和 KK, 的 数值 代入 (15) 式 ， 即 算出 质量 、 半 径 和 自转 角速度 在 5 个 演化 阶段 对 应 的 演化 时 标的 
数值 ， 计 算 结 果 见 表 1。 

R1 在 慢 引 力 收缩 阶段 金牛 T 星 的 质量 、 半 径 、 自 转角 速度 随 年 龄 的 演化 


Table 1 In the step of the slow gravitational contraction the evolution of mass, 


radius and rotational angular velocity of T Tauri with age 


m/mo R/R, w/w t (yr) 
0. 999 999 999 —0.279 479 5 12. 806 25 3. 250 2 x 10° 
0. 999 999 990 —0. 279 479 8 12. 806 23 3. 260 2 x 10° 
0. 999 999 900 —0. 279 482 2 12. 806 22 3. 261 0 x 10° 
0. 999 999 000 —0. 279 506 2 12. 800 17 3. 2609 x 10’ 
0. 999 990 000 —0. 279 746 87 12. 777 89 3. 261 1 x 10° 
质量 逐次 减少 半径 收缩 逐次 加 大 9 转角 速度 逐次 减 慢 时 标 
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(1) 恒星 诞生 后 开始 进入 快 引力 收缩 阶段 。 在 此 阶段 引力 收缩 占 主 导 地 位 ,质量 损失 可 以 不 计 ， 
自转 加 快 。 当 抛射 物质 时 进入 慢 引 力 收 缩 阶段 。 此 时 质量 损失 占 主导 地 位 ,引力 收缩 占 次 要 地 位 , A 
转角 速度 渐 慢 。 本 文 研 究 金牛 T 星 演化 快 引力 收缩 终止 ， 开 始 进 入 慢 引 力 收 缩 阶段 以 后 的 演化 进程 。 
表 1 给 出 的 数值 是 从 快 引力 收缩 开始 时 间 t=0 到 慢 引 力 收缩 开始 时 间 t= 3. 250 2 x 10 年 以 后 到 t= 
3.26 x 10 的 5 个 阶段 演化 的 数值 。 

(2) 质量 m 和 半径 R 联 立 演 化 方程 组 的 建立 是 在 表面 有 效 温度 7, 和 抛射 物质 速度 V 为 常量 的 情 
况 下 推出 的 。 表面 有 效 温度 视 为 常量 是 根据 第 1 节 第 (5) 条 ， 假 设 收缩 星 没 进 入 主 序 前 垂直 于 H-R 图 
的 横 坐 标 (光谱 型 坐标 上 的 有 效 温度 ) 的 路 径 而 演化 (Hayashi 演化 轨迹 ) 5 ， 但 这 只 是 近似 垂直 ， 只 
限于 没 进 入 主 序 前 期 的 演化 阶段 。 当 收缩 星 快要 进入 主 序 时 演化 路 线 向 左 弯 曲 而 不 垂直 于 H-R 图 的 
光谱 型 横 坐 标 ， 在 此 阶段 了 .为 常量 就 不 适用 。 此 外 ， 金 牛 T 星 抛射 物质 速度 了 视 为 不 变 。 在 慢 引 力 
收缩 阶段 ， 金 牛 T 星 的 质量 损失 主要 来 自 抛 射 物质 ， 因 能 源 只 靠 引力 收缩 能 而 尚 无 核能 ， 故 光子 辐 
射 (质量 辐射 ) 造 成 的 质量 损失 也 可 不 考虑 。 

(3) 表 中 给 出 的 数值 理论 上 是 合理 的 。 第 1 行 的 数值 表示 金牛 T 星 快 引力 收缩 阶段 终止 ， 进 入 慢 
引力 收缩 阶段 开始 时 间 z= 3. 25 x 10* 年 的 演化 质量 、 半 径 和 开始 减 慢 速度 的 数值 。 第 1 列表 示 因 抛射 物 
质 质量 逐年 减少 ; 第 2 列表 示 半 径 逐 年 收缩 尺度 加 大 ( 负 号 表示 半径 收缩 ); 第 3 列表 示 自 转角 速度 逐年 
减 慢 ; 第 4 列表 示 质 量 和 半径 的 演化 时 标 。 根 据 文 [12] 对 28 颗 金 牛 T 星 自转 的 观测 ， 自 转角 速度 的 数 
a 值 下 降 ， 但 下 降 较 少 。 本 文 表 中 自转 角速度 减速 的 变化 值 在 5 个 阶段 也 很 少 ， 这 说 明和 观测 一 致 。 表 
1 的 时 标 即 质量 从 0. 999 999 999m, 开 始 的 演化 时 标 也 是 半径 从 -0. 279 479 5R, 开 始 的 演化 时 标 。 


7 结 论 


本 文 研 究 的 金牛 了 星 是 其 它 类 型 的 金牛 类 星 的 典型 星 , 它 代表 大 多 数 金牛 类 型 星 的 普遍 性 质 。 在 
慢 引 力 收 缩 阶 段 ， 因 抛射 物质 使 质量 逐渐 减少 ， 引 力 收 缩 使 半径 逐渐 收缩 增 大 尺度 ， 自 转角 速度 逐渐 
减 慢 。 这 些 变化 是 理论 和 计算 的 结果 ， 有 待 观 测 证 实 。 但 目前 尚未 观测 到 金牛 T 星 的 自转 角速度 加 快 
的 事实 ， 正 如 前 节 讨 论 ， 根 据 文 [17] 对 28 颗 金 牛 T 星 自转 的 观测 ， 自 转角 速度 在 减 慢 。 
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Evolutional Time Scale of Mass and Radius of T Tauri and Its Influence on 
the Rotation under the Action of the Ejected Material and Gravitational 
Contraction in the Step of the Slow Gravitational Contraction 


Li Linsen 
(School of Physics, Northeast Normal University, Changchun 130024, China, Email; dbsd-lls@163.com) 


Abstract: The mass loss and radius contraction of T Tauri star under the role of the ejected material and 
gravitational contraction in the step of the slow gravitational contraction are studied. The evolutional equations 
determining the variation of mass and radius with time are established. The evolutional time scales of mass and 
radius are given. The influences of mass loss and radius contraction on the variation of the rotational angular 
velocity are given too. As example, the numerical effects of mass loss and radius contraction on the rotational 
relative angular velocity are estimated in Table. In addition the obtained results are discussed and concluded. 
Key words: T Tauri star; Slow gravitational contraction and ejected material; Evolution scales of mass and 


radius; Change of rotational angular velocity 


